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pL^ ' Abstract. We report on the June 2000 long (100 ks) BeppoSAX exposure that has unveiled above 10 keV a new 

very high energy component of the X-ray spectrum of rj Car extending to at least 50 keV. We find that the 2-150 
keV spectrum is best reproduced by a thermal + non-thermal model. The thermal component dominates the 
2-10 keV spectral range with kTh=5.5±0.3 keV and logNH/i=22.68±0.01. The spectrum displays a prominent 
iron emission line centred at 6.70 keV. Its equivalent width of 0.94 keV, if produced by the thermal source, gives 
^ ' a slightly sub-solar iron abundance ([Fe/H] =— 0.15±0.02). The high energy tail above 10 keV is best fitted by 

^\ I a power law with a photon index of 2.42±0.04. The integrated 13-150 keV luminosity of ~12 Lq is comparable 

OA ■ to that of the 2-10 keV thermal component (19 Lq). The present result can be explained, in the rj Car binary 

star scenario, by Comptonisation of low frequency radiation by high energy electrons, probably generated in the 
colliding wind shock front, or in instabilities in the wind of the S Dor primary star. It is possible that the high 
energy tail had largely weakened near the minimum of the 5.53 yr cycle. With respect to the thermal component, 
it probably has a longer recovering time like that of the highest excitation optical emission lines. Both features 
can be associated with the large absorption measured by BeppoSAX at phase 0.05. 
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CZJ 1. Introduction A similar behaviour was also found at other wavelength 

^ • r 1 bands, from radio to X-rays, which is commonly inter- 
im The peculiar southern obiect 77 Car is one or the most re- ri-ii ^ ■ t ■ 11 
-> ' .... >, , , . , preted m terms 01 a highly eccentric binary model com- 
markable variables m our Galaxy due to dramatic changes 1 r o t-> ^ 1 • • ^ 1 
. . Tir,-o. 1 11.1 posed 01 an S Uor-type very luminous primary star, and 
m its brightness. In 1843 it was the second brightest star 



o 
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^ • 1 , , rr 1 1 n 1- i i • i i ^n unsecu early type secondary star. The binary system 

S3 m the sky, then suuered a deep ladme down to the eighth . ^ , , ,,. . , 1-^1 1 1 

., , , ,ri^, f Tr-- interacts through colhdmg winds producing the observed 

■ - - ' magnitude by the end ot tlic 19th century (e.g. Viotti ^ i - v •■/ 

i,^r>r\ ^ 111- V Gars high temperature, luminous A-ray emission (e.g. 

1995). During the last century 77 Car was slowly and irreg- ^ggg^ Corcoran et al. 2001). 
ularly re-brightenmg up to the present V~5. Presently, 
according to the current distance estimates, rj Car has 

a bolometric magnitude around 5x10^ Lq (2x10^0 erg Recently, thanks to the BeppoSAX unique broad-band 

s-\ e.g. Hillier et al. 2002). A mass loss rate of lO'^/-^ x-ray coverage we were able to detect, for the first time, 

Mq yr^i or larger has been estimated from observations j] Car above 10 keV (Viotti et al. 1998; Viotti et al. 2002a, 

(e.g. Hillier et al. 2002, van Boekel et al. 2003, Pittard & Paper I). We reported the December 1996 observation 

Corcoran 2002, Andriesse et al. 1978). showing a 13-20 keV flux in excess with respect to the 

Optical spectroscopic observations unveiled a peculiar extrapolated 5 keV thermal spectrum that dominates the 

cyclic behaviour, showing regularly repeating excitation 2-10 keV range. The presence of a high energy tail was 

minima, with a period of 5.53 years (Damineh et al. 2000). confirmed by the following BeppoSAX observations of 31 

'7, Z a ~^ T' J -n u i -sr- ii- -i December 1999-2 January 2000 (Rebccchi et al. 2001). In 

bend offprint requests to: Roberto Viotti, e-mau: uvs- j \ i 

pace@rm iasf cnr it particular, in June 2000 a 100 ks exposure unveiled that 

* Based on space observations collected with the BeppoSAX t^e tail was probably non-thermal and extending to at 

X-Ray Astronomy Satellite which is a program of the Agenzia least 50 keV (Viotti et al. 2002b). In this work we anal- 

Spaziale Italiana with participation of the Netherlands Agency yse in detail the latter observation in order to investigate 

for Aerospace Programs. the origin of these very high energy photons, and compare 
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Fig. 1. The BeppoSAX MECS images of the 77 Car region on 2000 June 21-23, in three energy bands: (from left to 
right) 2-5.5 keV (softer), 6-7.5 keV (iron line), and 7.5-10 keV (harder). North is on the top and east to the left; 
the size of the image field is 54 arcmin. The second brightest source to the west of rj Car is the W-R star HD 93162 
(WR25). 



with the previous BeppoSAX observations of 77 Car. The 
results are summarised in Table 1. 

2. Observations 

The BeppoSAX satellite observed r] Car on June 21- 
23, 2000 with its Narrow Field Instruments (NFI). 
NFI include a Low Energy Concentrator Spectrometer 
(LEGS) (Parmar et al. 1997), and three Medium Energy 
Concentrator Spectrometers (MECS) (Boella et al. 1997) 
at the foci of four X-ray telescopes, a Phoswich 
Detector System (PDS) (Frontera et al. 1997), and an 
High Pressure Gas Scintillation Proportional Counter 
(HPGSPC) (Manzo et al. 1997). LEGS and MECS units 
have imaging capabilities and cover the 0.1-10 keV and 
1.5-10 keV energy ranges, respectively. PDS covers the 12- 
200 keV band, while HPGSP the 7-60 keV energy range. 
They are both coUimated instrument with a field of view 
of 1.3 and 1.1 degrees, respectively. The PDS collimators 
are rocked back and forth by about 3° to allow the simul- 
taneous monitoring of the source and background. During 
2000 HPGSP was switched off. rj Car has been observed 
for a nominal exposure time of 100 ks. The effective ex- 
posures were 80.8 ks for MECS no. 2 and 3 (MECS no.l 
was not in operation), and 33.4 ks for PDS. The corre- 
sponding count rates were, MECS: 0.639±0.003, and PDS: 
0.165±0.032 sec-i. The LEGS image will not be discussed 
here, for the paper is devoted to the analysis of the com- 
pact hard component which dominates beyond ^2 keV. 
The observation date corresponds to phase 1.457 of the 
spectral variation cycle of 5.53 y (Damineli et al. 2000), 
assuming $ = 1 for the 1998.0 minimum. 

Figure 1 shows the MEGS images of the region around 
rj Gar using photon events selected from three different 
energy bands, in order to trace the softer and harder X- 



ray sources in the field, and to map the regions emitting in 
the 6.7 keV iron line. It is evident in the figure that 77 Gar is 
the hardest and most luminous object in the region. Some 
other sources in the field present a non negligible residual 
in the 7.5-10 keV range and at the iron line, the most 
interesting one being the Wolf-Rayet star HD 93162/WR 
25 (WN6-t-04), to the West of rj Car. 

This observation, in particular, shows that 77 Gar is the 
main, and, most probably, unique contributor to the flux 
observed with the PDS instrument. As discussed in Paper 
I the PDS field of view also includes the X-ray pulsar IE 
10148.1-5937; but, according to Oosterbroek et al. (1998) 
and Tiengo et al. (2002) its spectrum is softer and the flux 
above 10 keV is much fainter than that of -q Car. Also the 
lack of intense diffuse emission in the high energy MECS 
band seems not to support the hypothesis of an extended 
very high energy emission due to the interaction of the 
winds from the massive stars in the Carina Nebula. 

3. Spectral analysis 

The MECS spectrum of 77 Car was extracted within a cir- 
cular region centred on the star with a radius of 4 ar- 
cmin. To take into account for possible contamination 
from the nebular X-ray emission, we extracted the back- 
ground from a coronal region around 77 Car with a 4.3 
arcmin inner radius and 11.3 arcmin outer radius. We 
also excluded from the selected area two circular areas 
of 4 arcmin in radius centred on the strong X-ray sources 
HD 93162 (WN6-I-04) and HD 93250 (03.5V). The back- 
ground subtracted spectrum of 77 Car and the local nebular 
spectrum normalised to the same area are shown in Fig. 2. 
The MEGS background subtracted spectrum of 77 Gar has 
been rebinned in order to have at least 20 source photons 
per energy bin. 
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Fig. 2. The BcppoSAX MECS background subtracted 
count rates (in keV~^) of r] Car in June 2000 com- 
pared with the local nebular background spectrum. 
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Fig. 3. Upper panel: Results of the fit of the June 2000 
BeppoSAX MECS background subtracted spectrum with 
a one temperature (5.51 keV) MEKAL model. The spec- 
tral model is extrapolated to the PDS energy range in 
order to show the high energy flux excess. The residuals 
are plotted in the lower panel. 



The background subtracted PDS data have been re- 
binned up to 50 keV following the standard procedure. 
Above 50 keV the data were binned in order to have a 2a 
count rate in each bin. As shown in Figure 3, r] Car has 
been detected at the >3a level up to 45 keV, and at the 2a 
level at higher energies. 77 Car had already been detected 
by PDS in December 1996 and December 1999, but only 
in the 13-20 keV range due to the shorter exposure times. 
A PDS upper limit was derived from the March 1998 ob- 
servation. Table 1 reports the 13-20 keV PDS count rates 
for the four BeppoSAX observations. 



The extracted MECS-hPDS spectrum of June 2000 was 
analysed using the XSPEC spectral fitting software pack- 
age. We started by analysing MECS alone in the 2-10 keV 
energy range which is dominated by the hard core compo- 
nent. The MECS spectrum was first fitted with a thermal 
model (MEKAL) leaving iron abundance free to vary. We 
obtained the following best fit values of the parameters: 
kTft=5.51±0.25 keV, NH/,=4.82±0. 12x10^2 cm'^, and a 
logarithmic iron abundance of [Fe/H]=-0.15±0.02. The 
reduced was 1.749 (for 64 degrees of freedom), partly 
due to the small instrumental energy mismatch at 6.7 keV 
discussed in Paper I. The source temperature is close to 
that measured in the MECS spectra of December 1996 and 
March 1998, and reported in Paper I. In June 2000 the 2- 
10 keV energy distribution is close to that of December 
1996, with nearly the same NH/j. As discussed in Paper I, 
NH/, was much larger in March 1998 during rj Car's egress 
from the X-ray minimum event ($ = 1.05; see also Figure 
6 below). 

As described in Paper I, we have also tried to 

model the observed MECS spectrum using an absorbed 
bremsstrahlung model and a Gaussian line. We obtained 
for the iron line a peak energy of 6.743±0.001 keV and 
an equivalent width of 0.964 keV, with values for kT/j 
(5.70±0.14) and NH,, ((4.91±0.10) x 10^^ cm^^) similar 
to those derived from the previous MEKAL model. As 
reported by Viotti et al. (2002b), The equivalent width 
of the iron line is consistent with that of December 1996 
(1.05 keV) and of December 1999 (1.01 keV), but smaller 
than that measured in March 1998 (see Table 1). The 
peak energy of the line, corrected for the -1-0.04 keV mis- 
calibration discussed in Paper I, suggests that the line 
is mostly emitted by the hot plasma, without any im- 
portant contribution from the 6.4 kcV fluorescent line. 
The reduced 1-96 in the range 2-10 keV. The 

integrated absorbed and unabsorbed fluxes in the 2-10 
keV energy range are: /2:=6.45x 10^^^ erg cm^^ s^^and 
/°=9.34x 10^^^ erg cm^^ s~^, respectively. The latter cor- 
responds to 19 Lq for 7] Car's distance of 2.6 kpc. 

Then, we have analysed MECS and PDS simultane- 
ously, allowing a PDS/MECS normalisation factor of 0.85, 
to take into account the miscalibration between the two 
instruments. This correcting factor was also used in Paper 
I, though not expHcitly indicated. As shown in Fig. 3, PDS 
counts are well above the extrapolation to higher energies 
of the MECS best fit thermal model. The difference be- 
tween observed and expected flux increases with energy, 
suggesting the presence of an additional very high energy 
component. 

In order to account for the high energy tail, we tried to 
fit the 2-150 keV spectrum with two thermal components. 
However, the fit of MECS-I-PDS with two freely varying 
temperatures always gives unrealistic results. Therefore, 
we fit MECS-I-PDS with frozen kT(l)=5.51 keV and 
NH(l)=4.82xl022 cm-2, and put constraints to kT(2) to 
be close to 10 keV. We thus derived: kT(2)=14±5 keV, 
and NH(2)=(543±22)xl022 cm'^, with a reduced of 
2.33 (Fig. 4). However, this result is not convincing, not 
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only because the fit runs below the observed flux above 30 
keV, but also because an unlikely very high NH is required 
for the higher temperature component. This suggests that 
the high energy spectrum of rj Car could be non-thermal. 

Therefore, we tried to reproduce the 2-150 keV spec- 
trum with a combination of thermal brcmsstrahlung (with 
frozen kT=5.70 keV) and power law components, having 
the same absorption column density (4.82x10^^ cm~^), 
and a Gaussian line at 6.7 kcV. The best fit model is shown 
in Figure 5. The non-thermal component has a best fit 
photon index of 2.42±.04 (reduced x^=2-41, for 79 dof). 
If produced in a shock, this photon index would imply a 
compression ratio for the shock front of x=3.1, thus the 
shock should be relatively weak (White & Chen 1994). The 
nonthermal component in the 13-150 keV range has an 
integrated energy flux of 5.51 xlO~^^ erg cm~^ (11.6 
Lq), a value comparable to that of the thermal 5 keV 
component. 

We have also attempted to fit the overall MECS+PDS 
spectrum with a single absorbed power law spectrum 
plus a Gaussian line. Unexpectedly, quite a good fit 
was obtained (excluding the photon events below 2 
keV) with a best fit photon index of 2.467±0.024, 
NH=(6.46±0.12)xl022, an Fe line centred at 6.733±0.007 
keV, and a reduced of 2.404. However, this solution ap- 
pears unrealistic, in particular because it would imply that 
the iron emission be mostly the 6.4 keV fluorescent line, in 
disagreement with the measured peak energy of the fea- 
ture (Ecor=6.70±0.01 keV). Therefore, the thermal plus 
power law model appears more likely to us. Incidentally, 
we notice that the two power law spectra have nearly the 
same slope, independently of the energy range considered. 

4. Discussion 

The high energy X-ray spectrum of rj Car is challenging, 
both because such an energetic and powerful source has 
so far never been seen in any other stellar-like source, and 
for no current model for non thermal emission seems at 
present applicable to the case of rj Car. 

In the T] Car's binary model it is assumed that the 
thermal (5 keV) X-ray emission is originated in the hot 
shocked gas produced by collision between the winds of the 
two stellar components (e.g. Pittard et al. 1998, Ishibashi 
et al. 1999). The shocked region should be placed at a 
distance from a few to many AU from the surface of the S 
Dor star during the highly eccentric orbital motion of the 
system. In this model, the plasma temperature is linked 
to the wind velocities, and would be unable to produce 
directly the highest energy photons seen in 77 Car. 

The power law spectrum exhibited by rj Car above 10 
keV indicates the presence of an additional non-thermal 
source. In principle, non-thermal emission can be ex- 
plained by several physical processes, most of which are, 
however, hardly compatible with current models of 77 Car. 
A plausible, but still difficult to accept process, would be 
inverse Compton scattering (IC) of low frequency photons 
by high energy electrons. In the case of r] Car a powerful 
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Fig. 4. Upper panel: Results of the fit of the BeppoSAX 
MECS and PDS (2-150 keV) spectrum of ry Car in June 

2000 with a two temperature (5.5 keV and 14 keV) 
MEKAL model. The residuals are plotted in the lower 
panel. 
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Fig. 5. The BeppoSAX MECS and PDS (2-150 keV) 

spectrum of rj Car fitted with a 5.70 keV thermal 
bremsstrahlung plus a power law spectrum with photon 
index 2.42, and a Gaussian line centred at 6.70 keV. 

photon source could be represented by the intense ultra- 
violet radiation from the stellar components of the binary 
system, probably from the hot secondary, since most likely 
the UV photons from the S Dor primary do not emerge 
from its dense wind. Stellar UV photons scattered by rel- 
ativistic electrons would finally emerge with energies in 
the X- and 7-ray range, and may carry away much of the 
energy pumped into the electrons at shocks. 

This process would require the presence of relativistic 
electrons with 7 ~ 100. According to White (1985) and 
Chen & White (1991) in the winds of early-type stars 
electrons can be accelerated to relativistic energies via the 
first-order Fermi mechanism in strong isothermal shocks 
generated by radiation driven instabilities in the winds. 
On the other hand, Jardine et al. (1996) have shown that 
in a colliding wind binary system, electrons can be acceler- 
ated up to relativistic energies in a current sheet, formed 
when the magnetic field carried in the stellar winds are 
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Table 1. Summary of the BeppoSAX observations during 1996-2000. Phase one corresponds to the spectroscopic 
minimum of 1998.0 
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forced together as the winds collide. In the case of r\ Car 
there is up to now no direct or indirect evidence of the 

presence of very high energy electrons. 

An indirect evidence could be non-thermal radio emis- 
sion, like that detected in galactic and extragalactic high 
energy X ray sources. However, non thermal radio emis- 
sion centered on t] Car's stellar core, has not been so far 
detected although we cannot exclude its presence. In fact, 
were it present it would lie below the surface that is opti- 
cally thick due to free-free emission. Actually, the whole 
central region of r\ Car is optically thick in radio wave- 
lengths to several arcscconds out. Since the non thermal 
emission should not extend that far away from the shock, 
it should be completely hidden by the thermal emission 
from the dense ionised stellar winds (Stephen White, priv. 
comm.). One should also consider that, as recently dis- 
cussed by Dougherty et al. (2003), in colliding wind binary 
systems the non thermal radio emission could be weak- 
ened by self-absorption, and inverse Compton cooling. 

Non-thermal radio emission has been detected in a 
number of massive binaries, suggesting the presence of 
relativistic electrons with power-law spectra, which could 
result in hard X- and 7-ray non-thermal flux via IC scat- 
tering (e.g. Benaglia & Romero 2003). Among them, the 
WR-t-0 binary system WR 140 bears some resemblances 
to the r\ Car system for its long period (7.9 yr), and high 
orbital eccentricity (e.g. Marchcnko et al. 2003), and for 
the presence of a variable IR excess attributed to dust for- 
mation near periastron (e.g. Williams et al. 1987). The X- 
ray emission of WR 140 is one order of magnitude smaller 
than that of r\ Car, but, like r\ Car, it suffers a strong 
absorption excess near phases 0.03 and 0.08 (Zhekov & 
Skinner 2000). If the physical process of production of 
non-thermal radio emission (that is of high energy elec- 
trons) in WR 140 is also working in t] Car, we would expect 
from 77 Car a very strong non-thermal radio emission. As 
discussed above, its absence is probably due to the heavy 
obscuration by the thermal absorption from the dense stel- 
lar winds. Indeed, the crucial point of the spectrum r\ Car's 
core radio emission remains at least at present unsolved. 

Self-Comptonisation is a mechanism proposed to ex- 
plain non-thermal emission in AGNs, as well as in galactic 
compact X-ray sources. As for instance discussed by Hua 
& Titarchuk (1995), quasi-power law X-ray spectra can 



be produced by diffusion of low frequency photons by op- 
tically thick relativistic plasma clouds. In this framework, 
we might assume that the spectrum of 77 Car above 10 keV 
be produced by scattering of the 5 keV thermal radiation 
itself by high energy (>10 keV) electrons present in the 
shocked region. This seems to us a more promising process 
as it does not require highly relativistic electrons as in the 
IC case. 

The study of the spectral variability along the 5.53 yr 
cycle could be a complementary way to tackle the problem 
of the origin of the high energy X-ray tail of r\ Car. Table 1 
reports the PDS coimt rates and the best fit parameters of 
the thermal component during the four BeppoSAX obser- 
vations. We recall that, according to the RossiXTE obser- 
vations, in December 1997 r] Car underwent a deep X ray 
eclipse, which is currently associated with the periastron 
passage of the suggested binary system (Ishibashi et al. 
1999). The X ray eclipse lasted until the end of February 
1998. Hence, the BeppoSAX observation of mid March 
1998 was made during j] Car's egress from the eclipse. 

Fig. 6 compares the BeppoSAX MECS spectra ob- 
served during the four epochs. For the sake of homogene- 
ity with the other observations, only MECS units 2 and 
3 have been considered for the December 1996 spectrum. 
It is evident in the figure that beyond the iron feature the 
four spectra nearly overlap, also suggesting that the highly 
absorbed March 1998 spectrum had about the same tem- 
perature and unabsorbed flux as in the other three epochs. 

As for PDS, it turns out that on three cases the PDS 
count rate in the 13-20 keV range was the same within 
the data uncertainty. As shown in Table 1, in March 1998 
77 Car was not detected with PDS. The 3(t upper limit of 
0.15 s~^ is equal to the flux measured in the other three 
observations, but we think that this coincidence is acci- 
dental, and that in March 1998 the high energy tail was 
weak. Most probably, this component followed the trend of 
the thermal one, and largely weakened during the 1998.0 
X-ray eclipse. Our March 1998 observation suggests that, 
while the thermal component had already recovered its 
luminosity, although with a residual large absorption, the 
high energy tail was still weak. It should have recovered its 
luminosity at a later time (but well before phase 1.37). It 
may be noticed that a long recovering time was also found 
in the visual spectrum for the high excitation emission 
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Fig. 6. Plot of the BeppoSAX MECS23 spectra of 77 Car 
during December 1996- June 2000. The phases of the 5.53 
yr cycle are indicated. Note the larger flux of the iron line 
in March 1998 ($ — 1.05) when the spectrum was strongly 
absorbed. At energies >7 keV the four spectra basically 
overlap each other. 

lines (e.g. [Ne iii]), which, as discussed in Paper I, were 
still very weak in May 1998. The limited time-coverage 
does not allow us to conlude whether there is a physical 
link between the two phenomena, but it would be worth- 
while investigating it in the future. 

5. Conclusion 

We have presented the first in-depth analysis of the spec- 
trum of the 77 Car system above 10 keV, based on a long 
exposure BeppoSAX observation. The power law best fit 
suggests a non-thermal origin of the hard X-ray tail. The 
integrated 13-150 keV luminosity (^12 Lq), is compara- 
ble to the luminosity of the 2-10 keV thermal component 
(19 Lq), and suggests the presence of a very effective for- 
mation process. 

Two models have been considered in the light of the 
proposed binary nature of 77 Car. In one model the high 
energy tail is produced by inverse Compton scattering of 
the UV stellar photons by relativistic electrons produced 
in the wind of the primary star, or in the shocked col- 
liding wind region. Alternatively, the high energy pho- 
tons are produced by self-Comptonisation of the thermal 
5 keV emission from relativistic electrons with energies 
much lower than in the previous case. Though, as sug- 
gested by the referee, the weakness of the high energy tail 
in March 1998 could be better explained as high inverse 
Compton cooling of the reletivistc electrons during peri- 



astron passage, when the colliding winds shock is closer to 
the stars. Neither model has so far enough support from 
observations, like detection of non-thermal radio emission 
from the central source, a point which would deserve new 
very high resolution radio observations. 

A crucial point for understanding the nature of the 
high energy tail of 77 Car, would also be to investigate the 
slope of the spectrum at higher energies and to measure 
the high energy cutoff of the spectrum, which is related 
to the energy of the scattering particles. Our PDS obser- 
vations allow us to determine a lower limit to the cutoff 
energy of ^50 keV. INTEGRAL observations might allow 
measurement of the X-ray spectrum of 77 Car above 100 
keV, and to determine up to what energy the power law 
spectrum is extending. 

Finally, new high energy (>10 keV) X-ray observa- 
tions, e.g. with the foreseen ASTRO-E satellite, of ry Car 
near the periastron passage of the system will provide 
a clue of where the non-thermal source is located. Our 
March 1998 upper limit might suggest a recovering time 
slower than that of the thermal source, as observed in 
the optical spectra in the high energy emission lines. It 
would important to investigate the physical link beteen 
the two phenomena, and whether that behaviour is asso- 
ciated with the high NH^ value found also at the eclipse 
egress. 
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